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I. INTRODUCTION 

In order to describe the motion of massive bodies in post- 
Newtonian approximations to General Relativity and other metric 
theories of gravity, one is led to define certain relativistic 
potentials. They simply add to the Newtonian potential, 
inducing perturbations on the Keplerian motion of a body in the 
gravity field of another body. Well-known examples of such 
perturbations are the Einstein and Lense -Thirring rotations of 
the orbit in inertial space. 

The question is whether the same relativistic potentials 
also induce perturbations in the rotational motion of a nearly 
rigid body and, if they do, how can these perturbations be 
observed? 

Recently [l] the author has developed an approach for an 
approximate description of the motion of extended bodies. The 
approach is based on the local hydrodynamical equations of 
motion for a perfect fluid in the PPN approximation [2] and, 
as such, is well adapted to the description of perfect fluid, 
weakly gravitating bodies, a class to which a large fraction of 
systems of astronomical interest belong. 

After a review of the formalism (Section II) , the main 
result is obtained in Section III: given a body A with spin 

A 

vector S and traceless part of the tensor of inertia (or 

A 

quadrupole moment) I a g in the field of a body B with mass m 0 , 
the potential which depends on the orientation of the ellipsoid 



S 


* 





of inertia of A has the structure 


V = - R*S + f 


jA jA jA 

■f- R c R B + ^ jSS- V 6 * h J8. R a Rg . 


( 1 . 1 ) 


where fi is defined exactly in (3.14); f, g, and h are scalar 
functions; R is the position vector of the Newtonian center- 

A 

of -mass of A with respect to the one of B; R is the unit 


A* A 

vector; and R = dR/dt. 

Thus the modification of the classical torque-inducing 
potential : 


3 G m B _A 
2 R 3 aB 


R R q 

a 8 


( 1 . 2 ) 


basically the second term in (1.1), consists of a potential 
inducing a relativistic spin-precession about the direction of 
ft, and of two new potentials which induce velocity-dependent 
torques. 

In the case of the Stanford Gyro Relativity Experiment 
[3.1, it is known that the detection of the first few terms in 
ft and of the classical gravity-gradient term induced by (1.2) 
is within the capability of the readout system. It is natural 
to ask whether the velocity-dependent torques could also be 
measured in the experiment. The answer is negative; the reason 
is that both g and h depend linearly on the product of velocity 
of body A (the gyroscope) and the velocity of body B (the Earth), 
measured in any inertial frame. Because of the large mass of 


2 





the Earth in comparison to that of the gyro, we can always 
choose an inertial frame in which the Earth is at rest, so 
that g = h = 0. 

This analysis then, in particular, confirms that the only 
gravitational causes which may affect the rotational motion of 
the gyroscope are the first few terms of the relativistic spin 
precession and the Newtonian gravity gradient (1.2). 

Throughout this report use will be made of the notion of 
a nearly rigid body. While in an exact relativistic context 
Born's concept of rigid rotation has no meaning and dynamical 
conditions (see [41) must instead be used, in the approximate 
framework used here the Born kinematical condition can never- 
theless be assumed valid in an approximate sense. Also, 
deviations from Born's rigid rotation will be assumed as small 
as terms of second post-Newtonian order. 

Units such that G = C = 1 are used throughout. Greek 
indices range from 1 to 3, Latin indices from 0 to 3. 
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II. FORMALISM 


In this section we proceed to summari7e some of the 
definitions and the main results in References 1 and 5. 

In post-Newtonian hydrodynamics for a perfect fluid it 
is convenient to introduce the following four independent 
variables and an equation of state as the basic set of quanti- 
ties in a given coordinate system: 
p mass density 

v^ coordinate three-velocity (2.1) 

f = II (p) equation of state relating the internal 
specific energy density I to p. 

In terms of the basic set (2.1), the following quantities 
are now introduced: 

components of the metric tensor: 


g ab “ g ab (p ' V ” ' 


( 2 . 2 ) 


gravitational red-shift factor: 


u °= - <- g 0 0 - %6 v a v 8 ' 2g 0a v a rlS ' 


(2.3) 


coordinate four-velocity of the matter-stream: 


v • = (1, v a ) , 


(2.4) 


pressure defined by the equation of state and the 
restriction to adiabatic processes: 


p : = pdlf/dlnp , 


(2.5) 


i>v' 

1 


j 
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energy-momentum tensor of a perfect fluid: 


T ab : = p |_(u 0 ) 2 (1 + H + p/p) v a v b + (p/p)g a j . (2.6) 

The energy momentum tensor (2.6) is entirely defined in 
terms of the basic set (2.1), provided that one knows how to 
explicitly specify — in the given coordinate system — the metric 
tensor in the form (2.2), or, in other words, how to solve the 
field equations. For a large number of systems of astrophysical 
interest a "slow motion" approximate representation of the 
type (2.2) is possible in a neighborhood of the material source 
and in harmonic coordinates by means of an iterative solution 
method of a reduced Einstein equation [1,6,7] or by means of 
the method originally devised by Chandrasekhar [8l or Fock [9]. 

In the following the metric, and correspondingly (2.3) 
and (2.6), will be considered known in the form (2.2) in the 
first post-Newtonian approximation. 

Besides the Newtonian potential at the point (x, t) 


U: * ; c=const e'lx-xT 1 d 3 *' • 


(2.7) 


one defines the post-Newtonian potential 


W i = f. . p'v' Ix-x'l A d^x' , 

a t=const K a 1 1 ’ 


( 2 . 8 ) 


which is a gravitational analog of the vector potential in 
classical electrodynamics. 

Let 6 and y be PPN parameters (6 = y = 1 in General 
Relativity) , and 
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(2.9) 


= p/^g u 


P(1 + j 


+ 3yU) 


the "conserved energy density" [21. 

Given an isolated (in the sense of [10]) two-body system 
(A, B) , 1 at any point inside the spatial support of A, say, one 
can uniquely decompose each potential into the sum of an 
"internal" and "external" part by simply restricting the domain 
of integration in (2.7) and (2.8) to the spatial support o^ 
body A (internal potential) and of body B (external potential) . 
Let denote the external potentials on A, and E^, E , 

the nonrelativistic kinetic, potential, and internal energy 
densities, respectively, of a given body, as measured by a 
comoving observer. 

Define for each body an inertial mass-density 


P 


I 



+ E 
£ 



( 2 . 10 ) 


along with the corresponding inertial mass 


m : = /, , p_d x , 

body 


( 2 . 11 ) 


and a gravitational mass-density 


P G : > (jll + lE p /P r ) 


(n : - 4 p - y -3) 


( 2 . 12 ) 


along with the corresponding gravitational mass 


■” (G) = = 'body P G ^ ' 


(2. 13) 


1. The same considerations can be readily extended to any 
iso] a ted system containing N > 2 bodies. 
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which in some metric theories of gravity may be different from 
the inertial mass (2.11) [11], 

Let the bodies (A,B) be homogeneous and in rigid rotation, 
a, b the coordinates of the respective Newtonian centers of mass, 
and the relative positions in A, B with respect to a,b, 

respectively. 

If one defines the PPN linear momentum of A 
( * : ~ d/dt; #a : = 9/3x a ) 

p a : * i 1 + 7 ) V« ' 2(1 + Y) U + 


+ 7 V ; B p ' V e |X ' X ' ' >a6 d3x d3x ' + 


+ (1 + 2y) / A pv^ U^ A) d 3 x , 


(2.14) 


and the PPN force on A, 


- [i + + ^ ^ h 


B P^ix-x' I , a d 3 x d 3 x ' + 


. a d3 - 


-1 


+ (1 + 2y) f A p J"g p ' (a-£ A + t.|^) |x-x' | , a d 3 x ^x' - 


- 2(1 + Y ) / A cv v W (A) v>cl d 3 x + 


+ 7 ; A pv 6 J'b p '''y | x ' X ' 1 *a6y d3x d3x '- 


(2.15) 


2. In the terms of post-Newtonian order the densities of 
gravitational and inertial mass can be replaced by p because 
the difference is negligibly small in this approximation. 
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then one can show [ 1] that the local equations of motion for 
the material source 

T ab ; b = 0 (2.1*' 

imply 



modulo quantities of second post-Newtonian order. 

Expanding the gradient of the potentials in (2.15) into 
multipoles and considering only the monopole term, one can 
construct out of (2.17) a Laqrangian which proves to be identical 
to the so-called "Einstein-Infeld-Hof fmann" Lagrangian [12,13], 
valid for spherically symmetric field-singularities, provided 
one identifies the mass (2.11) or (2.13) of a perfect fluid, 
rigidly rotating body with the "mass" of a field-singularity. 

Considering one further term in the multipole expansion of 
the potentials introduces a relativistic dependence of the 
Lagrangian on the spin of the bodies. Use of the theory of 
perturbation of the initial data (i.e., the orbital elements 
of the osculating Keplerian ellipse) leads to a secular precession 
of the direction of the orbital angular momentum or, equivalently, 
of the longituc e of the ascending node of the relative orbit. 

Since ihe spin-dependent Lagrangian has been i ained from 
the equation of the translational motion (2.17), it cannot be 
used in order to study the rotational motion of the bodies of 
the system. Caporali and Spyrou [5] have, however, shown that 


4 >’ 
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if one introduces a generalized definition of the spin which 
includes terms of post-Newtonian order, the equations of motion 
of the spin can be obtained, for each body, from the theorem of 
conservation of the total angular momentum of the system and 
the equation of motion of the orbital angular momentum, which 
is already known from orbital analysis. In this way they were 
able to show that the resulting equation of motion of the spin 
of each body was the same as that which could be obtained by 
standard techniques from the Lagrangian of the translational 
motion. The proof can in principle be repeated at any order of 
the multipole expansion. It can then be stated that at any order 
of the multipole expansion in (2.17) the corresponding Lagrangian 
contains for each body the correct dependence on the three 
coordinates of its center of mass and on its three Euler angles. 

A Lagrangian with the same dependence on .hese six degrees 
of freedom per body was obtained by Barker and O'Connell using a 
quantum approach to the relativistic equations of motion of 
extended bodies, assumed to be described by the energy-momentum 
tensor of spin 1/2 particles [14]. 

The following section is dedicated to the study of the 
relativistic rotational motion of the bodies in the isolated 
system. The strategy will be as follows: in the multipole 

expansion of the relativistic potentials, the terms depending 
on the Euler angles of a given body will be singled out from those 
depending on the translational degrees of freedom to obtain the 
potentials entering the rotational (Euler) equations of motion of 
that body. 
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III. THE POTENTIAL OF THE ROTATIONAL MOTION 

In this section I ta,’ 1 sketch the calculation which leads 
us to recognize the dependence of the relativistic Lagrangian 
corresponding to (2.17) on the Euler angles of a given body A, 
say (the formulae valid for B can be obtained from those of A 
by simply interchanging the labels) . Once this dependence 
has been established, the Euler equations give straightforwardly 
the time-evolution of the spin of the body. 

Let a, b be the coordinates of the Newtonian center-of-mass 
of A, B , respectively. Replacing in the force (2.15) the quantity 
|x - x'|~\ by its Taylor expansion in a neighborhood of (x = a, 
x' = b) , one obtains a collection of terms which can be classi- 
fied according to whether they are "torque -free" or "torque- 
inducing" . 

Let us first examine the torque-free terms; that is, those 

A 

corresponding to a potential of the form ft • S , for some ft, 

A 

S being the spin of A, defined below. Let 

S* £ - 5. 

be the relative position-vector in A. 

In (2.15) one recognizes that the terms which have nonzero 
spin dependence are 

(1 + 2y) a p f B p' | x-x' j -1 a d 3 x d 3 x' (3.1) 

and 

-2(1 + y) ; A pv g f B p* v 0 ' | x— x 1 | — 1 d 3 x d 3 x' . (3.2) 


10 



»*■ - • - 1 ' - 1 • •*• ' *..*.Vi4*i ftijfll 
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The latter contributes through 


-2(1 + y) b*j^p j^p'lx-x’l" 1 Q dx d 


3 ^3 , 

v * x* 


(3.3) 


and 


-2(1 + r) / A P ? A - / B P‘ £’ B Ix-x'l" 1 >a d 3 x d 3 x' 


(3.4) 


Define the spin according to classical mechanics 


“ e * o8 r i p <-1 «s d3x 


I = A, B. 


and the tensors as 


_B r j-B .B .3 

I v....v 4 p ^V....^v d x ' 

1 n 1 n 


n B , ZB r B r B ,3 

P av 1 '"v ' f B p ^vi . . . d x 
1 n 1 n 


(n > 2) 


( 3 . 5 ) 


which are symmetric in the indices v i‘** v n * 

The spin-dependent part of (3.1) or (3.3) is the gradient 
of 

D -1 

t-xj jB 3 m+x; R 

D3 n! v x ... v "• " u 

K n=2 n 


1 . C A , R « (-l) n t B 

7 paY S Y B _ 3 y n! I v 1 ...v 


v x.-V 


) (3.6) 


contracted with a, b, respectively. 

Only the traceless part of the I's, defined by Reference 


15, 3 


Vj . . . v. 


int (n/2) 

= . E i-» k 

n k=2 


iysi.6 . 6 i b 

2 n (Vj^ 2 ... v zj c _^ v 2 j c v 2k+i • • * v n * ’ * * p n P n 


(Sj 


(3.7) 


5. int (-) denotes the integer part of parentheses 
2 2 

on indices denote complete symmetrization . 
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is involved in (3.6). 


Also, introducing the operator 


rB ^ _(-l) n jB 9 (n) 


* - £ 


n _2 Vj....v n Vj . . . . v n 


and the notation 


<a I B , VR _1 > : = 


e _ a i B (3 .r” 1 ) , 
jiaB a y 


(3.8) 


(3.9) 


the spin-dependent part of (3.1) is the gradient of 


i±2i («. 6JL4 - <a I B , VR->». S A 
O o n 3 


(3.10) 


and likewise for (3.3). 

The spin-dependent part of (3.4) is given by the gradient of 


-2(1+y) [j < S B R-> 


I e paB 


s* £ 

n=2 


(- 1 ) 

n! 


civ. 



(m-1) 

v,.. 



(3.11) 

Again, only the traceless part P of the P’s in (v^...v ), 
defined according to the same rule as in (3.7) , contributes in 

(3.11). Defining the vector-valued operator p as in (3.8) 


P B , = £ p B 3 t») 

n=2 ~ 1 n i n 


(3.12) 


and using the same notation as (3.9), with replacing al, the 
spin-dependent part of (3.6) can finally be written as the 
gradient of 

- < 1+1r >[} S e 3 aS R '' - </> 

Introduce the total and reduced mass 


vv » ' ■ 

' ' ’ * **» -fst* ‘ I 
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f«*?Srv 


> ^ 4H 


» ro, + m B 
O A B 


y - m A m B /m 0 


and center-of-mass coordinates: 


m_a ■ nR ■ -m_b^ . 
A a a B a 


The torque-free potential is finally given exactly defining 

m. 


( 3 . 13 ) 


- + 1 + U ^ [S B - 3 (S B • R) R] 

- — iL < R i B , VR _1 > + (1+y) < ^ R i B J> B , VR _1 >| 

The torque-free potential is then V_-, « - 0 • S A . In an 

Tr — •*'* 

inertial frame locally comoving with A, the torque-free 
potential induces a precession of the spin of A about the 
direction of il: 

dS A - o y c A 

art* ii £ • 

The first two terms in the right-hand side of (3.14) have 
already been computed, though by means of different techniques 
( 16,171. The lowest order contributions to £_ from the last two 
terms can be written defining 


(3.14) 


< - »A 


m B R 


(3.15) 


One obtains: 




m_ 


B 


m. 


+ 1+ y) (R * R - 


4 Q* R) - ^4 S B - 3(S B -R)R 

‘ 2R 3 L ~ J 


+ R“ 2 0[(f./R) 3 1 . 


(3.16) 
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r > 


* 


C 


This expression agrees with a formula obtained by Barker and 
O'Connell when m << m [18] . 

Let us consider now the torque-inducing potential; that is, 
the one depending on the instantaneous orientation of the 
ellipsoid of inertia of A w. r.t. the position of B and vice versa. 
An exact formula in the same fashion as (3.14) must represent 
the coupling of all the multipoles of A with those of B, and it 
is not very illuminating, in general. In the following I shall 
consider for the bodies the multipoles up to and including the 
quadrupole, neglecting in the potential terms of order 0[(£/R) 3 ]. 

It is convenient to define 



m 


<(G)aB 


B 


(G) 


+ i B m 
(G)aB A 


(G) 



(3.17) 


where the label G indicates that the mass and the quadruple 
moment correspond to the gravitational mass density (2.12). 

The contribution to the torque-inducing potential coming 
from the first term in (2.15) is readily seen to be 

[l + ( * 2 + ^ 2) ] M aB 8 aB R_1 + R ~ 2 °t ( S/ R)3 ] 

modulo terms of second post-Newtonian order. 

Likewise, the contribution from the second term in (2.15) is 


(i ’ 2?) 4 M aB 3 a8 R " X + R " 2 °n?/R) 3 l , 

while the third integral gives no contribution, at this order. 


14 



Finally, the contributions from the fourth and fifth 


integral in (2.15) are, respectively, 

-U+y> i*bM ag 3 2 & R” 1 + R“ 2 0[ (S/R) 3 ] 


and 


"wv 3 «8uv r + R ' 2 °f <5/ R > 3 1 • 


Writing 


M &v > * v K 


! = ■ R " ' 


Sv K v 


the torque- inducing potential is: 


v t '[ 1 + il r (i2 + * 2 > + t - < 1+ *> k± * i k a V 6 ae - 56 A»] 


. - -fab,(6 ( 


_ x - 6 W + R ~ 2 °[ K/» 3 l • 

R R 


(3.18) 


A^ain, making use of the center-of-mass coordinates, we 
note that 


• * A A* 

R a = R R a + R R a 


so that 


*a6 Ve R *RR + R *RR 


(3.19) 


and 


*oB R a k B = ** l Rf, + 2 R R + r2 • 


(3.20) 
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The first term in the right-hand side of (3.19) and (3.20) 
contributes to the first term in the right-hand side of (3.18). 
By contrast, igg and Igg introduce in the torque- inducing 


potential a nonclassical dependence of the orientation of the 
ellipsoi 1 of inertia of one body with respect to the instantaneous 
position and velocity of the other body. 

If the body A has a mass much smaller than that of B, the 
last term in (3.18) tends to zero and can be neglected. In fact, 
it will always be possible to choose an inertial system in which 
b = 0. Accordingly, the velocity-dependei _, torque-inducing 
potentials will not be involved in the equation of the .-otational 
motion of A. 
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IV. CONCLUSION 
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We have proved that in the post-Newtonian approximation the 
potential governing the rotational motion of an extended body 
of mass m A , spin S A and quadrupole moment in the relativistic 

f p\ 

field of a mass m' 1 has the structure 

O 

ro (C) 

V = - R.S A + -L- l'°> R a R v (« ve - 6 R v R 6 > M/n» o ♦ 


+ [1 + F (6 , y» orbital parameters)] 3li G ^ m^ G ^/2 R^ . (4.1) 

RR B 


The potential (4.1) is then the post-Newtonian generaliza- 
tion of (1.2). 

An exact expression for (2 has been derived in Section III 
in terms of the differential operators (3.8), (3.12). 

The second term in (4.1) is responsible for a nonclassical 
dependence of the orientation of the ellipsoid of inertia of 
one body with respect to the instantaneous position and velocity 
of the other. When the body A has a much smaller mass than that 
of B, the effect of this term on the rotational motion of the 
body tends to zero. 

The fact that the nodal line of an orbiting body, defined 

by the intersection of the equatorial plane and orbital plane, 

rotates on the orbital plane is confirmed by the third term in 

(4.1). Only the frequency of rotation needs to be redefined. 

In practice, this consists of multiplying the classical frequency 

by the sum of 1 and the numerically small function, of order 
2 

(V/C) contained in the parentheses. 
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It is worth noting that the framework developed in 
Reference 1 leads to the association with each body of not 
only a gravitational mass but also of a countable set of 
infinitely many gravitational multipole momenta. This fact 
is a consequence of the existence of a gravitational mass 
density (2.12) which may differ from the inertial one (2.10) 
in some metric theories of gravitation. 

Observable consequences on the relativistic rotational 
motion of a small, nearly rigid body can be induced only by 
the torque-free term -ft*S under the form of a precession of 
the spin S about the direction of ft. The feasibility studies 
of the Stanford Gyro Relativity Experiment (see, e.g., [3]) 
have pointed out that this precession can be observed as a 
variation in the orientation of the London magnetic moment of 
a spinning quartz ball in geodetic motion about the Earth. 
Careful evaluation of the attainable level of accuracy has 
shown the feasibility of a measure of the leading part (3.16) 
of the vector ft. 
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